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Abstract. We consider the Schrodinger operator H = (i V + A) 2 in the space 
L2(R 3 ) with a magnetic potential A created by an infinite straight current. 
We perform a spectral analysis of the operator H almost explicitly. In par- 
ticular, we show that the operator H is absolutely continuous, its spectrum 
has infinite multiplicity and coincides with the positive half-axis. Then we 
find the large-time behavior of solutions exp(—iHt)f of the time dependent 
Schrodinger equation. Equations of classical mechanics are also integrated. 
Our main observation is that both quantum and classical particles have al- 
ways a preferable (depending on its charge) direction of propagation along the 
current and both of them are confined in the plane orthogonal to the current. 



1. Introduction 

They are very few examples of explicit solutions of the Schrodinger equation 
with a magnetic potential. Probably the only ones are a constant magnetic field, 
B(x,y,z) = B , (sec, e.g., [3]) and, in the two dimensional case, a magnetic field 
localized at the origin, B(x,y) = B S(x,y) where S(x,y) is the Dirac function (see 
[1]). The solution is expressed in terms of Hermite functions in the first case and in 
terms of Bessel functions in the second case. Here we suggest a third example of an 
explicitly solvable Schrodinger equation. Actually, we consider the magnetic field 
B(x,y,z) created by an infinite straight current. Physically, this case is opposite 
to the case of an (infinitely) thin straight solenoid considered in [1] by Aharonov 
and Bohm where the field is concentrated inside the solenoid. 

Suppose that the current coincides with the axis z and that the axes x, y and z 
are positively oriented. According to the Bio-Savart-Laplace law (see, e.g., [5]) 

B(x,y,z) = a(-r~ 2 y,r- 2 x,0), r = (x 2 +y 2 ) 1/2 , 

where \a\ is proportional to the current strength and a > (a < 0) if the current 
streams in the positive (negative) direction. The magnetic potential is defined by 
the equation 

B(x, y, z) — curl A(x, y, z) 

and can be chosen as 

(1.1) A(x,y,z) = -a(0,0,lnr). 

Thus, the corresponding Schrodinger operator in the space Z/2(R 3 ) has the form 
H = H y = -dl -flg + (id, -7lnr) 2 , 7 = ea, 
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where e is the charge of a quantum particle of the mass m = 1/2 and the speed of 
the light c = 1. 

Since magnetic potential (1.1) grows as r — > oo, the Hamiltonian Zf does not fit 
to the well elaborated framework of spectral and scattering theory. Nevertheless, we 
perform in Section 2 its spectral analysis almost explicitly. To be more precise, we 
reduce the problem to an ordinary differential equation with the potential 7 2 In 2 r 
(let us call it logarithmic oscillator). We show that the operator H is absolutely 
continuous, its spectrum has infinite multiplicity and coincides with the positive 
half-axis. Then we find in Section 3 the large-time behavior of solutions exp(—iHt)f 
of the time dependent Schrodinger equation. In Section 4 we integrate equations of 
classical mechanics. Our main observation is that positively (negatively) charged 
quantum and classical particles always move in the direction of the current (in the 
opposite direction) and are localized in the orthogonal plane. 

A detailed presentation of the results of this note can be found in [6] . 

2. Spectral analysis of the operator H 
Let us first consider a more general magnetic potential 

(2.1) A(x,y,z) = (0,0, A(x,y)) 

with an arbitrary (we disregard here domain questions) real function A(x, y) which 
tends to infinity (cither +oo or — oo) as r = (x 2 + y 2 ) 1 ^ 2 — ► oo. The corresponding 
Schrodinger operator is 

H = -A + (id z +eA(x,y)) 2 , 

where A is always the Laplacian in the variables (x,y). Since A does not depend 
on z, we make the Fourier transform $ = <fr z in the variable z. Then the operator 
H = <f>i7$* acts in the space L 2 (M. 2 x R) as 

(Hu)(x,y,p) = {h(p)u)(x,y,p}, 

where 

(2.2) h(p) = -A + (p-eA(x iy )) 2 . 

Here pGR (the momentum in the direction of the z-axis) is the variable dual to z 
and the operator h(jp) acts in the space L 2 (M. 2 ). 

Since A(x,y) — > oo or A{x,y) — ► — oo as r — > oo, the spectrum of each operator 
h(p) is positive and discrete. Let A„(p), n = 1,2,..., be its eigenvalues, numerated 
in such a way that A„(p) are analytic functions of p. The spectrum of the operator 
H, and hence of H, consists of the sets (branches) covered by the functions X n (p), 
n = 1,2,..., as p runs from — oo to oo. This is similar both to the cases of 
the constant field where A„(p) = \eBv\(2n + 1) + p 2 and to the periodic problem 
where the role of the momentum p in the direction of the z-axis is played by the 
quasimomentum (see, e.g. [4]). Thus, the general Floquct theory implies that 
the spectrum of H is absolutely continuous, up eventually to some eigenvalues of 
infinite multiplicity. Such eigenvalues appear if at least one of the functions A n (p) 
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is a constant on some interval. Then this function is a constant for all p € K. On 
the other hand, if, say, the function — eA(x, y) is semiboundcd from below, then 

(2.3) lim inf (p - eA(x, y)) 2 = oo 

as p — > oo, and hence linip^oo X n (p) — oo for all n. Thus, we have the following 
simple result. 

Theorem 2.1. Suppose that A(x,y) is a semi-bounded function which tends either 
to +00 or to —00 as r — > 00. Then the operator H is absolutely continuous. 

Note that the Thomas arguments (see, e.g., [4]) relying on the study of the 
operator-function h(p) for complex p are not necessary here. 

The problem may be further simplified if A(x, y) — A(r). Then we can separate 
variables in the polar coordinates (r, 9) . Denote by TL m the space of functions 
u(r)e lm0 where u e L 2 (R;rdr) and m = 0, ±1,±2, ... is the orbital quantum 
number. Then 

00 

(2.4) £ 2 (K 2 )= n m . 

m— — oo 

Every subspace H m is invariant with respect to the operator h(p). The spectra 
of their restrictions h m (p) on TL m consist of positive simple eigenvalues \ m ,i{p) < 
Kn,2(p) < • • • , which are analytic functions of p. We denote by ip m .i(r,p), tp m ,2{r,p) 7 ■ ■ ■ 
the corresponding eigenfunctions which are supposed to be normalized and real. 

Let us return to the operator H with the potential A(r) = — alnr. In this case 
operator (2.2 ) equals 

h(p) = -A + ln 2 (e p r 7 ). 



Since H 7 u = H- 7 u, it suffices to consider the case 7 > 0. It is convenient to transfer 
the dependence on the momentum p into the kinetic energy and to introduce the 
parameter a — e p / 7 G (0, 00) instead of p. Let us set 

(2.5) K(a) = -a 2 A + 7 2 ln 2 r, 

and let w(a), (w(a)f(x, y) = af(ax, ay), be the unitary operator of dilations in the 
space L 2 (M 2 ). Then 

(2.6) w*(a)h(p)w(a) = K(a), a = e ph . 

We denote by /i m ,n(a) and (j) 7 n,n{r, a) eigenvalues and eigenfunctions of the re- 
strictions of the operators K(a) on the subspaces Tt m - It follows from (2.6) that 
Vm,n( a ) = \n,n{p) and <Pm,n( a ) = w * (^V'm.n (p) ■ Actually, decomposition (2.4) 
is needed only to avoid crossings between different eigenvalues of the operators 
h(p). It allows us to use formulas of perturbation theory (see, e.g., [4]) for simple 
eigenvalues. We fix m and omit it from the notation. 

The next assertion is quite elementary but plays the crucial role in the following. 

Lemma 2.2. For every n, we have that p' n (a) > for all a > 0. 
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Indeed, analytic perturbation theory shows that 

(2.7) f j,' n (a) = {K' n (a)(p n {a),4> n {a)) = 2a \V<j) n (x,y,a)\ 2 dxdy. 

Jr 2 

This expression is obviously positive since otherwise <j) n (x, y, a) — const. 

The next assertion realizes an obvious idea that the spectrum of K(a) converges 
as a — > (in the quasiclassical limit) to that of the multiplication operator by 
7 2 In 2 r, which is continuous and starts from zero. 

Lemma 2.3. For every n, we have that lim a _^ Mn( a ) = 0- 

Since the function lnr is not semibounded, relation (2.3) is not true in our case. 
Nevertheless, taking into account the kinetic energy, we obtain the following result. 

Lemma 2.4. For every n, we have that lim^oo n n {a) = oo. 

In terms of eigenvalues X n (p) of the operators h(p), Lemmas 2.2 - 2.4 mean that 
Ki(p) > f° r all p £ IR and linip^-oo \ n (p) = 0, limp^oo A„(p) = oo (for 7 > 0). 

Let A n be multiplication operator by the function A ra (p) in the space L2O&). It 
follows from the results on the function A„ (p) that the spectrum of A„ is absolutely 
continuous, simple and coincides with the positive half axis. Let us introduce a 
unitary mapping 

00 

* : L 2 (R+ x R; rdrdp) L 2 (R) 

n=l 

by the formula (*/)„(p) = J °° f(r,p)ip n (r,p)rdr. Then 

OC 

(2.8) *$ff$***=0A n 

n=l 

(of course Jf = iJ m and A„ = A n m ) , and we obtain the following 

Theorem 2.5. The spectra of all operators H m and H are absolutely continuous, 
have infinite multiplicity and coincide with the positive half axis. 

As a by-product of our considerations, we have constructed a complete set of 
cigenfunctions of the operator H. They are parametrized by the orbital quantum 
number m, the momentum p in the direction of the z-axis and the number n of an 
eigenvalue \ m ,n{p) of the operator h m (p) defined by formula (2.2) on the subspace 
H m - Thus, if we set 

«m,n, P (r, -z, 0) = e iF e"" 9 i,»(r,p), 
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3. Time evolution 

Explicit formulas obtained in the previous section allow us to find the asymp- 
totics for large t of solutions u(t) = exp(— ii/i)uo of the time dependent Schrodinger 
equation. On every subspace with a fixed orbital quantum number m, the problem 
reduces to the asymptotics of the function u(t) = cxp(—iH m t)u . Below we fix m 
and suppose that 7 > 0. 

Assume that 

(3.1) ($u )(r,p)=ip n (r,P)f(p), 
where / e Cft°(R). Then it follows from formula (2.8) that 

/OO 
-OO 

The stationary points of this integral are determined by the equation 

(3.3) z = X' n (p)t. 

Since X' n (p) > 0, the equation (3.3) has a solution only if zt > 0. We need the 
following information on the eigenvalues /i n (a) of the operator (2.5). 

Lemma 3.1. For every n, we have that lim a _> a^! n {a) = 0. 

Indeed, it follows from equation (2.7) that a/i' n (a) < 2/i n (a). Therefore it remains 
to use Lemma 2.3. 

Lemma 3.1 means that linip^-oo A^(p) = 0. The following conjecture is physi- 
cally quite plausible and is used mainly to formulate Theorem 3.3 below in a simpler 
form. 

Conjecture 3.2. For every n, we have that A" (p) > for all p € K and lim^oo (p) = 
00. 

Therefore equation X' n (p) = v has a unique solution p n = ip n (v) for every v > 0. 
Clearly, 

(3.4) \UMvW n (v) = l. 

Let 3>„(v) = (p n (v)a — A rl (< f o„(u)), 9(v) = 1 for v > 0, 9{v) = for v < and 
±i = e ±7rl//2 . Applying to the integral (3.2) the stationary phase method and 
taking into account identity (3.4), we find that 

(3.5) 

u(r, z, t) = e^'M^r, ^{z/t^z/t) 1 ' 2 /(^(zA))^)- 1 / 2 ^^/*) + Uoo (r, z, t), 
where 

(3.6) lim ||«oo(-,t)|| = 0. 

Note that the norm in the space L2 (K+ x R) of the first term in the right-hand side 
of (3.5) equals ||uo||- The asymptotics (3.5) extends of course to all / e L 2 (R) and 
to linear combinations of functions (3.1) over different n. Thus, we have proven 
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Theorem 3.3. Assume that Conjecture 3.2 is fulfilled. Suppose that 7 > 0. Let 
u(t) = ex.p(—iH m t)uo where u n satisfies (3.1). Then the asymptotics as t — > ±00 
of this function is given by relations (3.5), (3.6). Moreover, if f 6 Cq°(M.) and 
=Fz > 7 then the function u(r,z,t) tends to zero faster than any power of \t\~ x as 
t — > ±00. 

Conversely, for any g G L 2 (R+) define the function uq by the equation 



Let us consider the motion of a classical particle of mass m = 1/2 and charge e in 
a magnetic field created by potential (2.1) where A(x, y) = A{r), r = (x 2 + y 2 ) 1 ^ 2 . 
We suppose that A(r) is an arbitrary C 2 -function such that A(r) = o(r _1 ) as r — > 
and |*4(r)| — > 00 as r — > 00. The solution given below is, to a large extent, similar to 
the Kepler solution of equations of motion for a particle in a spherically symmetric 
electric field. However, in the electric case the motion is always restricted to a plane, 
whereas in the magnetic case it is confined in the plane z = but the propagation 
of a particle in the ^-direction has a non-trivial character. We proceed here from 
the Hamiltonian formulation. An approach based on the Newton equations can be 
found in [6]. 

Let r be a position of a particle and p be its momentum. Let us write down the 
Hamiltonian 



in the cylindrical coordinates (r,tp,z). In the case (2.1) where A{x,y) = A(r), we 



where p rj p v and p z are momenta conjugated to the coordinates r, ip and z. Since 
H(r,p) does not depend on (p and z, the momenta p v {t) and p z (t) are conserved, 
i.e., Pipit) = M (M is the moment of momentum with respect to the z-axis) and 
p z (t) = P (the magnetic momentum in the z-direction). Therefore Hamiltonian 
equations read as 




4. Classical mechanics 



H(r,p) = (p 2 -eA(r)) 2 



have that 



H(r,p) = (p 2 r + r- 2 pl) + (p z -eA(r)) 2 , 



(4.1) 



r'(t) 




where 



(4.2) 



V(r) = M 2 r- 2 + (P- eA(r)) 2 



and 



(4.3) 



if'(t) = 2Mr(t) 



-2 



(4.4) 



z'(t) = 2(P-eA(r(t))). 
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It suffices to solve the system (4.1) since, given r(t), the solutions of equations 
(4.3) and (4.4) are constructed by the formulas 

(4.5) ip(t) = ip(0) + 2M [ r(s)- 2 ds 

Jo 

and 

(4.6) z(t) - «(0) = 2 / (P - eA(r{s)))ds. 

Jo 

The solution of the system (4.1) is quite similar to the solution of the Kepler problem 
although in our case the effective potential energy (4.2) depends additionally on the 
momentum P in the z-dircction. In the solutions of the quantum problems, it is 
reflected by the fact that, for electric spherically symmetric potentials, the variables 
can be separated (in the spherical coordinates), whereas in our case the operators 
h(jp) depend on p. 

Thus, to solve (4.1), we remark that 

(4.7) A^r'itf + V(r{t)) = K, 

where K = 4~V(0) 2 + M 2 r(0)~ 2 + 4~ 1 z'(0) 2 is a constant kinetic energy of a 
particle. Clearly, (4.7) is the equation of one-dimensional motion (see, e.g., [2]) 
with the effective potential energy V(r) and the total energy K. It admits the 
separation of variables and can be integrated by the formula 

(4.8) t = ±4 J{k- V(r)j dr. 

Note that V(r) — > oo as r — > and r — > oo. Let r m ; n and r max be the roots of 
the equation V(r) = K (r m ; n and r max are the nearest to r(0) roots such that 
Tnim < r(0) < r max ). It follows from (4.8) that the function r(t) is periodic with 
period 

/•''max, ._l/2 

(4.9) T = 8 (K-V(r)J dr 

J r m i„ 

and r m i n < r(t) < r max . One can imagine, for example, that on the period the 
function r(t) increases monotonically from r min to r max and then decreases from 
r max to r m ; n . Thus, we have integrated the system (4.1) and (4.3), (4.4). 

Theorem 4.1. In the variable r a classical particle moves periodically according 
to equation (4.8) with period (4.9). The angular variable is determined by equation 
(4.5) so that ip(t) is a monotone function of t and ip(t) — (fat + 0(1), where ipo = 
2MT~ 1 f Q r(s)~ 2 ds, as \t\ —> oo. The variable z(t) is determined by equation (4.6). 

According to equation (4.4) a particle can move in the direction of the current 
as well as in the opposite direction. Nevertheless one can give simple sufficient 
conditions for the inequality 



(4.10) 



± (z(t + T) - z{t)) > 
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(for all t). Indeed, it follows from the Newton equation r"(t) = — 2V'(r(i)) (which 
is a consequence of (4.1)) and expression (4.2) that 

r"(t) = 4M 2 r- 3 (t) + 4eA'(r(t))(P - eA(r(t))). 

Using equation (4.4), we see that 

2ez'(t) = (r"(t) - 4MV 3 (t))i'(r(t))- 1 . 

Integrating this equation and taking into account periodicity of the function r(t), 
we see that, for all t, 

2e(z(t + T) - z{t)) = f T r"(s)A'(r(s))- 1 ds - AM 2 C r(s)~ 3 A'(r(s))~ 1 ds 
Jo Jo 

(4.11) = [ r'(t) 2 A'(r(t))- 2 A"(r(t))dt-4M 2 [ r(«)- i U / (r(*)r 1 d*- 
Jo Jo 

Let us formulate the results obtained. 

Theorem 4.2. The increment of the variable z on every period is determined by 
equation (4.11). In particular, if ±eA'(r) < and ±eA"(r) > for all r, then 
inequality (4.10) holds. In this case z(t) = z$t + 0(1) with zq = T^ 1 (z{T) — z(0)), 
±z > 0, as \t\ — > oo. 

In particular, for potentials A(r) = — a\nr and A(r) = —ar a where a £ (0,1), 
inequality (4.10) holds if ±ea > 0. Note that in these cases the fields B(x,y,z) = 
A'(r)r^ 1 (y, —x, 0)tend to as r — > oo. 

It follows from equation (4.4) that if, say, eA'(r) < and the point r cr is 
determined by the equation p = eA(r cr ), then z(t) increases for r(t) £ (r cr ,r max ) 
and decreases for r(t) £ (r m i n , r cr ). Of course, it is possible that r cr < r m ; n . In this 
case, z(t) always increases. Let us discuss this phenomena in more details on our 
leading example A(r) = — alnr. Then r cr = e~ p l~< where 7 = ea. The points r m ; n 
and r max are determined from the equation 

V(r) = M 2 r- 2 + ln 2 (e p r 7 ) = K. 

The function z(t) is increasing for all t if r cr < r m j n or, cquivalently, V(r cr ) > K 
and r cr < r(0). The first of these conditions is equivalent to M 2 e 2P ^ > K or, 
since in view of (4.4) e 2P h = r (0)- 2 e z ' ( ^^ , to 

M 2 r(0r V'(°^ > 4-V(0) 2 + M 2 r(0)" 2 + 4- 1 z'(0) 2 . 

Thus, z'(Q) should be a sufficiently large positive number (z'(0) < is definitely 
excluded). In this case the condition r cr < r(0) which is equivalent to z'(0) > 
is automatically satisfied. Note finally that always r cr < r max , that is the function 
z(t) cannot be everywhere decreasing (this is of course also a consequence of Theo- 
rem 4.2). Indeed, inequality r cr > r max is equivalent to V(r cr ) > K and r cr > r(0). 
The first of them require that z'(0) > while the second require that z'(0) < 0. 

Thus, positively (negatively) charged classical and quantum particles always 
move asymptotically in the direction of the current (in the opposite direction). In 
the plane orthogonal to the direction of the current classical and quantum particles 
are essentially localized. 
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